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I .  INTRODUCTION 

The  work  described  In  this  report  was  conducted  between  l  October  1975 
and  31  March  1977  under  Defense  Advanced  Research  Projects  Agency  contract 
DAAB  07-76-C-1342,  M.I.  Bell,  principal  Investigator,  F.  Rothwarf,  contract 
monitor.  This  Investigation  continued  the  efforts  begun  under  a  previous 
contract  (DAAB  07-74-C-0470)  to  elucidate  the  physical  mechanisms  responsible 
for  the  pyroelectric  effect  in  ferroelectrtcs ,  to  discover  such  limitations 
as  may  exist  on  the  performance  of  these  materials  in  Infrared  Imaging 
(vldicon)  systems,  and  to  obtain  theoretical  guidance  in  the  search  for 

materials  which  would  yield  Improved  detector  performance.  These  objectives 

1-3 

have  been  met  through  the  development  of  a  generalized  molecular  field 
theory  of  ferroelectrlcity  (GMFT)  which  has  provided  new  insight  Into  the 
nature  of  ferroelectric  phase  transitions  and  has  led  to  a  better  under¬ 
standing  of  the  observed  limitations  on  the  performance  of  the  materials 

4 

currently  used  as  vldicon  targets  .  A  brief  account  of  the  basic  elements 
of  this  theory  has  been  prepared  for  publication  and  Is  contained  in 
Section  II  of  this  report. 

An  underlying  assumption  of  the  GMFT  Is  the  existence  of  elementary 
dipoles  which  persist  (although  In  a  disordered  state)  at  temperatures 
well  above  that  of  the  ferroelectrlc-to-paraelectrlc  transition,  even 
in  so-called  displacement-type  ferroelectrtcs.  In  order  to  determine 
whether  effects  attributable  to  the  presence  of  such  dipoles  could  be 
observed  In  the  Raman  scattering  spectra  of  displacement-type  ferroelectrtcs, 
a  study  of  KNbO^  was  begun  In  collaboration  with  A.M.  Qulttet  and  M.  Lambert"* 
and  completed  with  the  help  of  M.  Krauzman.  The  results  of  these  experiments, 
which  provide  excellent  support  for  the  assumptions  on  which  the  GMFT  Is 
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based,  were  presented  In  a  previous  report,^  and  a  revised  version  of  that 

manuscript  has  since  been  published  In  the  Physical  Review^.  A  reprint 

of  this  work  Is  provided  In  Appendix  A. 

Further,  experimental  verification  of  the  GMFT  has  been  obtained 

from  measurements  of  the  electric  field  and  temperature  dependence  of  the 

dielectric  constants  of  several  ferroelectrlcs .  From  these  measurements 

one  can  obtain  the  temperature  dependence  of  the  coefficients  which  appear 

8  12 

In  the  Devonshire  expansion  of  the  free  energy  .  The  GMFT  predicts  * 

that  these  coefficients  should  depend  linearly  on  temperature,  and  our 

experimental  results,  presented  in  Section  III  of  this  report,  confirm 

this  prediction.  An  account  of  these  experiments  will  be  presented  at  the 

March  1978  meeting  of  the  American  Physical  Society,  and  an  abstract  of 

that  talk  Is  provided  In  Appendix  B. 

Our  assessment, based  on  the  GMFT,  of  the  potential  performance  of 

4 

ferroelectrlcs  as  vldlcon  targets  led  us  to  conclude  that  materials 

significantly  better  than  those  now  In  use  are  not  likely  to  be  found 

within  the  well-known  families  of  ferroelectrlcs  (e.g.  perovskltes  or 

trlglyclne  sulfate  derivatives)  or  by  minor  modifications  of  these 

materials.  We  therefore  turned  our  attention  to  situations  In  which 

the  arguments  of  the  GMFT  would  not  apply  and  one  might  hope  to  avoid 

the  limitations  Imposed  by  the  theory.  One  such  case  Is  that  of  the 

9 

Improper  ferroelectrlcs  (e.g.  the  rare-earth  molybdates),  where  an 
analysis  based  on  phenomenological  models  for  the  free  energy  shows 
that  the  performance  of  such  a  material  as  a  vldlcon  target  depends  on 
temperature  In  a  significantly  different  way  than  the  performance  of 


_ 
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a  proper  ferroelectric.  The  calculated  performance  (slgnal-to-t»lse 
ratio)  of  a  proper  ferroelectric  target  decreases  rapidly  as  the 
Curie  temperature  is  approached,  but  when  an  Improper  ferroelectric 
is  used  the  performance  is  expected  to  approach  a  maximum  at  the 
Curie  temperature.  These  results  were  presented  to  a  meeting  of  the 
Materials  Research  Council  of  the  Defense  Advanced  Research  Projects 
Agency  in  July,  1976  (see  Appendix  C  for  an  abstract  of  that  report) . 
Section  IV  of  this  report  contains  the  details  of  these  calculations 
together  with  preliminary  experimental  results  which  both  demonstrate 
the  theoretically  predicted  behavior  and  indicate  that  the  increase  in 
performance  exhibited  by  an  improper  ferroelectric  near  its  transition 
temperature  can  be  great  enough  (more  than  a  factor  of  15)  to  make  it 
competitive  with  the  best  of  the  materials  now  in  use  as  vidlcon  targets. 

Finally,  Section  V  examines  several  potential  applications  of 
pyroelectrics  in  the  areas  of  energy  conversion  and  storage  which, 
while  not  strictly  within  the  scope  of  the  present  investigation,  are 
of  considerable  current  interest  and  depend  for  their  success  on  the 
development  of  materials  having  many  of  the  same  properties  as  are 
required  in  thermal  imaging  systems.  For  each  application  a  figure  of 
merit  is  developed  which  permits  an  assessment  of  the  relative  merit  of 
the  various  ferroelectric  materials  which  might  be  used  in  that  application. 
Results  from  the  GMFT,  from  the  phenomenological  theory,  and  from  the 
experimental  literature  are  used  to  evaluate  and  compare  these  figures 
of  merit,  and  several  promising  directions  are  identified  for  future 
studies  aimed  at  achieving  efficient  thermal  to  electrical  energy 
conversion  using  ferroelectrics . 
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The  most  recent  and  fruitful  theory  of  ferroelectric  phase  transitions 
is  based  on  the  concept  of  the  lattice  dynamical  instability,  or  soft  mode, 

introduced^  by  Cochran  and  Anderson.  Older  theories  employing  the 

2 

molecular-field  approximation  describe  the  ferroelectric  phase  transition 

as  an  order-disorder  transformation  involving  the  ordering  of  an  array  of 

permanent  dipoles,  as  in  the  Weiss  theory  of  ferromagnetism.  Despite  a 

3 

recent  renewal  of  interest  in  this  approach  ,  it  would  appear  not  to  be 
applicable  to  a  large  class  of  ferroelectrics,  including  the  perovskites, 
where  the  structure  determined  for  the  paraelectric  phase  has  each  atom 
located  at  an  inversion  center,  thus  precluding  the  existence  of  permanent 
dipoles.  In  addition,  the  phase  transitions  predicted  by  the  conventional 
molecular-field  theory  (MFT)  are  thermodynamically  of  second  order,  while 
many  observed  transitions  are  first  order.  The  purpose  of  this  communication 
is  to  show  that  neither  of  these  objections  should  be  regarded  as  ruling 
out  the  possibility  of  an  adequate  description  of  ferroelectricity  based  on 
a  molecular-field  theory. 

4 

With  regard  to  the  first  objection,  it  has  been  pointed  out  that  the 

atomic  sites  may  be  inversion  centers  only  in  some  average  sense.  For 

example,  studies^  of  the  scattering  of  electrons  and  x>~rays  by  BaTiO^  and 

4 

KNbO^  reveal  patterns  of  diffuse  scattering  which  can  be  interpreted  as 

arising  from  partially  correlated  displacements  of  the  metal  ions  from 

the  centers  of  the  oxygen  octahedra.  The  observed  scattering  is  accounted 

for  by  assuming  a  relatively  simple  pattern  of  displacements  resulting  in 

the  formation  of  chains  of  correlated  displacements  approximately  10  to  25 

4 

unit  cells  in  length  ,  Only  when  averaged  over  distances  large  compared  to 
the  chain  length  does  the  crystal  have  the  macroscopic  structure  determined 
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by  x-ray  diffraction.  More  recently,  similar  diffuse  scattering 

has  been  seen  in  neutron  scattering  experiments^  on  BafTiO^,  KTaO^, 

KTa  Nb  0_  (KTN) ,  and  SrTiO-.  Additional  evidence  for  disorder  in  the 
x  1 "X  3  ^ 

perovskites  nas  been  obtained  from  Raman  scattering  measurements^  and 

Q 

from  studies  of  optical  second-harmonic  generation  . 

9 

Although  NMR  measurements  have  shown  that  the  displacement 

4 

correlations  are  dynamic  rather  than  static  as  originally  proposed  ,  this 
is  not  an  obstacle  to  the  use  of  a  molecular-field  model.  As  noted  by 
Yamada  et  al^,  the  unit-cell  dipole  moment  which  appears  in  such 
calculations  can  be  interpreted  as  an  instantaneous  value,  corresponding 
to  one  of  several  sets  of  equilibrium  atomic  positions  arranged  symmetrically 
about  the  average  structure.  The  free  energy  obtained  by  time  averaging 
would  then  have  the  same  form  as  if  the  dipoles  were  static. 

If  we  take  seriously  the  possibility  that  the  paraelectric  phases 
of  materials  which  undergo  ferroelectric  phase  transitions  may  be  polar 
in  the  sense  of  containing  permanent  electric  dipole  moments,  we  must 
re-examine  the  order-disorder  transitions  predicted  by  the  MFT.  For 
the  sake  of  this  discussion,  3nd  without  any  real  loss  of  generality, 
we  will  consider  an  array  of  N  elementary  dipoles  of  moment  jj,  ,  free  to 
assume  any  orientation  in  space,  which  interact  via  a  molecular  (local) 
field  E^.  The  usual  statistical  argument  then  gives  for  the  macroscopic 
polarization  P  at  temperature  T 

P(T)  -  NyL(x),  (l) 

where  L(x)  *  cothx  -  l/x  and  x  ■  fiE^/kT  (k  is  Boltzmann's  constant). 

For  small  values  of  x,  Eq.  (1)  can  be  solved  to  give  the  local  field 
E^  as  a  power  series  in  P. 
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It  has  been  customary  to  assume  that  E ^  has  the  Lorentz  form 

=  E  +  XP  (2) 

where  E  Is  the  applied  (macroscopic)  electric  field,  and  the  Lorentz 
factor  X  has  the  value  4tt/3  for  cubic  or  isotropic  media.  One  can  then 
solve  for  E  =  9G/SP  and  obtain  the  free  energy, 

G  =  %A0(T  -  Tn)  P2  +  &qT  P4  +  J  C0T  P6+...,  (3) 

which  has  the  Landau-Devonshire  form^  (i.e.  a  power  series  in  the 
order  parameter  P) ,  with 

Ac  =  3k/N/i2 

B  =  3 A.  /5N 2 1?  (4) 

o  o  r 

C  =  99A  /175N4  u 
oo 

and 

T  =  X/A  .  (5) 

oo 

In  this  context, the  most  important  fact  about  the  free  energy  (3)  is 

L 

that  B,  the  coefficient  of  the  P  term,  is  positive  for  all  finite 

temperatures,  so  that  the  MFT  always  predicts  a  second-order  phase 

transition^.  Yet,  as  will  now  show,  the  polar  nature  of  the  para- 

electric  phase  offers  a  means  of  avoiding  this  difficulty. 

12 

It  has  long  been  known  that  the  Lorentz  local  field  (2)  is  not 

suitable  for  the  description  of  liquids  composed  of  polar  molecules. 

With  this  local  field  the  MFT  predicts  a  ferroelectric  ordering 

12 

in  many  commm  liquids  at  room  temperature.  Water,  for  example  , 
would  order  at  Tq  =  1140  K.  Among  the  numerous  attempts  to  avoid  this 
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difficulty  one  finds  a  particularly  interesting  and  elegant  approach 
13 

due  to  Onsager  .  He  showed  that  in  a  polar  medium  the  field  which 
acts  to  align  an  elementary  dipole  with  the  applied  field  is  not  given 
by  Eq.  (2)  but  by 

*1  "  *  +  2&  <6> 

If  e  >  >  1,  Eq*  (6)  gives  E^^E  +  fcx^P,  and  since  the  inverse 

“12  2 

susceptibility  is  given  by  X  *  -  9  G/9P  ,  an  expansion  of  the  form  (4) 
for  G  requires  that  contain  terms  involving  every  odd  power  of  P: 

- E  +  \,p  +  ^ip3  +  v5+-'-  (7) 

The  Onsager  local  field  (6)  takes  into  account  only  the  dipole-dipole 

14 

coupling  and  ignores  any  short-range  interactions  .  For  this  reason 
we  treat  the  appearing  in  Eq.  (7)  as  adjustable  (temperature- 
independent)  parameters,  rather  than  use  the  specific  (temperature- 
dependent)  values  predicted  by  the  Onsager  theory  and  given  implicitly 
by  Eq.  (6).  The  generalized  local  field  (7)  leads  to  a  free  energy 


G  -  %Aq(T  -  Tq)  P2  +  %Bq(T  -  Tt)  P4  +  |  Cq(T  -  Tj)  P6  +.  . . ,  (8) 

where  the  coefficients  Ao  ,  BQ ,  CQ  , . . .  are  given  by  Eq.  (4),  and  the 
characteristic  temperatures  by 


T  =■  X  /a 
o  o  o 


Ti  -  H/Bo 


T2  "  /Co  * 


(9) 


and  so  forth.  The  most  striking  feature  of  this  generalized  molecular 
field  theory  (GMFT)  is  that  for  T^  >  Tfl  the  second  Devonshire  coefficient, 

(10) 


B  -  B0 (T  -  Tt)  , 
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Is  negative  in  the  vicinity  of  T0,  and  one  obtains  a  first-order  phase 
transition. 

In  comparing  the  predictions  of  Eq .  (8)  with  experimental  results, 

we  find  that  in  the  two  cases  where  the  temperature  dependence  of 

B  has  been  measured,  the  behavior  predicted  by  Eq.  (10)  is  actually 

observed.  Table  I  lists  the  values  of  the  coefficients  \  ,  B  and  the 

o  o 

temperatures  T0,  obtained  in  BaTiO^^  and  triglycine  sulfate  (TGS) 

Note  that  in  BaTiO-  we  find  T,  >  T  ,  as  required  for  a  first-order 
J  l  o 

transition,  while  we  have  T,  <  T  in  TGS  where  the  transition  is 

i  o 

second  order.  This  experimental  data  enables  us  to  make  two  quantitative 
checks  of  the  GMFT.  First,  using  Eq .  (4),  we  find  for  the  maximum 
(saturation)  polarization 

Nn  -  (3Aq/5Bo)*,  (ID 

2 

and  the  data  of  Table  I  yield  values  of  17  and  4.9  |jC/cm  for  BafTiO^ 

and  TGS,  respectively.  The  result  for  BaTiO^  is  somewhat  smaller  than 
2  15 

the  26  |jC/cra  found  experimentally  ,  but  that  for  TGS  is  in  excellent 

17  2 

agreement  with  the  experimental  result  of  4.3  |£/cm  .  Second,  Eq.  (4) 
also  gives  the  relation 

N  =■  A2/5kB  (12) 

o  o 

for  the  number  of  independent  dipoles  per  unit  volume.  We  find  densities 
21  -3  20  -3 

of  1.99  x  10  cm  and  4.6  x  10  cm  for  TGS  and  BaTiO^,  respectively. 
Comparing  these  results  to  the  unit  cell  volumes,  we  obtain  0.8  and  33  unit 
cells  per  elementary  dipole,  indicating  that  while  each  unit  cell  in  TGS 


behaves  as  a  statistically  Independent  dipole,  the  effective  dipoles 

in  BaTiO^  must  consist  of  larger  regions  in  which  the  polarization  is 

\  4 

highly  correlated.  The  linear  correlations  proposed  by  Comes  et  al 

to  account  for  the  diffuse  x-ray  scattering  have  a  length  of  10  -  25 

unit  cells  and  would  satisfy  this  prediction  of  the  GMFT  quite  well. 

A  major  virtue  of  the  present  model  is  its  ability  to  describe 

first  order  transitions,  and  so  special  attention  should  be  paid  to 

those  quantities  which  are  peculiar  to  such  transitions,  i.e.  the 

difference  between  the  Curie-Weiss  temperature  Tq  and  the  actual 

transition  temperature  T  ,  the  discontinuity  P  in  •pentaneous  polarization, 

c  sc 

and  the  latent  heat  of  the  transition  W.  Unfortunately,  accurate  calculation 

of  these  quantities  from  the  model  is  frustrated  by  a  lack  of  experimental 

information  concerning  the  coefficients  X^,...  which  determine  the 

higher-order  nonlinearities  of  the  local  field.  The  fact  that  the 

discontinuity  in  spontaneous  polarization  at  T  in  BafTiO.  is  707.  of 

c  3 

the  saturation  polarization*5  indicates  that  these  nonlinearities  are 
significant  even  at  the  transition  temperature.  We  have  found  that  this 
saturation  behavior  can  be  taken  into  account  in  an  ad  hoc  but  effective 
way  by  replacing  the  power  series  (7)  by  Its  lowest-order  nontrivial 
Pade  approxtmant. 


X  P 

o 

l-(Xx/Xo)P2 


(13) 


This  expression  reproduces  the  first  two  terras  of  (7)  for  small  P  and 

also  insures  that  if  X^  >  0  the  polarization  will  saturate,  since  then 
2 

P  <  X./X  for  all  finite  values  of  the  local  field. 

I  o 


11 


Using  Eq.  (13)  for  the  local  field,  we  find  for  the  spontaneous 
polarization  at  the  transisiton  Pgc/NM,  =0.6,  in  reasonable  agree¬ 
ment  with  the  experimental  result  cited  above.  We  also  obtain 

Tc  -  To  =  5K,  somewhat  less  than  the  experimental  value^5  of  8K. 

1 8 

The  latent  heat  (integrated  over  the  interval  Tq,  Tc)  is  cal¬ 
culated  to  be  21  cal /mol,  while  the  measured  value  is  47  cal /mol. 
Clearly,  the  predicted  transition  is  not  as  strongly  first  order  as 
that  observed  experimentally,  but  the  agreement  is  still  quite 
satisfactory.  Further  comparison  of  this  theory  with  experiment  and 
a  detailed  discussion  of  the  local  field  (13)  will  be  presented 
elsewhere . 
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I.  Experimental  values  of  parameters  appearing  in  the  free 

energy  given  by  Eq .  (8) . 


BaT103  (Ref.  15)  TGS  (Ref.  16) 


5  Cgs) 

7.6 

389 

13  Cgs) 

0.18 

110 

378 

322.4 

448 

247 
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ABSTRACT 


The  dielectric  constants  of  the  paraelectric  phases  of  barium 
titanate  and  triglycine  sulfate  are  measured  as  functions  of  tempera¬ 
ture  and  applied  electric  field,  and  the  results  used  to  determine 
the  temperature  dependence  of  the  coefficients  appearing  in  the 
Devonshire  expansion  of  the  free  energy.  Each  of  the  coefficients 
measured  is  found  to  have  a  linear  temperature  dependence,  a  result 
which  is  in  qualitative  agreement  with  a  recently  proposed  generalized 
molecular  field  theory  of  ferroelectricity.  Certain  anomalies  appear 
in  the  data  for  TGS ,  however,  and  these  will  require  further  study. 
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I.  INTRODUCTION 


According  to  Devonshire's  thermodynamic  theory  of  ferroelectrlcity  , 

if  a  ferroelectric  crystal  has  a  center  of  inversion  above  the  Curie 

temperature,  then  its  free  energy  function  G,  at  zero  stress,  can  be 

2 

expanded  in  terms  of  P  .  Assuming  that  in  the  ferroelectric  phase  the 
spontaneous  polarization  lies  along  the  z-direction,  and  that  electric 
fields  are  applied  only  along  this  direction  so  that  P=Pz>  Px=Py=0,  we 
get 


G=>^AP  +  k BP 


♦K 


(i) 


where  the  coefficients  A,3,C . are  functions  of  the  temperature. 

The  applied  electric  field  is  given  by 

E=  =  AP  +  BP3  +  CP5  + .  (2) 

and  the  dielectric  stiffness  (reciprocal  susceptibility)  is 

2 

X~l  =  =  A  +  3BP2  +  5CP4  + -  (3) 

3P* 

In  order  to  account  for  the  Curie-Wetss  law  behavior  of  the  dielectric 
constant,  A  Is  taken  to  be  a  linear  function  of  temperature  A*A0(T-T J . 
All  the  higher-order  Devonshire  coefficients  are  assumed  to  be  tempera¬ 
ture-independent  in  most  discussions  of  the  thermodynamic  theory. 

2 

In  1955,  Drougard  et  al  studied  the  polarization  dependence  of 
X  ^  in  the  paraelectric  phase  of  BaTTlO^  and  found  that  B  is  a  linear 
function  of  temperature  B=Bo(T-T1).  More  recently.  Bell  and  Raccah3,4 
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proposed  a  generalized  molecular  field  theory  of  ferroelectriclty  (GMFT) 
which  predicts  that  all  the  Devonshire  coefficients  should  be 

linear  functions  of  temperature.  Consequently  the  present  study  was 
undertaken  to  Investigate  experimentally  the  temperature  dependence  of 
the  Devonshire  coefficients.  We  performed  our  experiments  on  barium  tltanate 
(BaTlO^)  and  trlglyclne  sulfate  (TGS)  ,  since  the  phase  transition  Is 
generally  believed  to  be  first  order  In  the  former  and  second  order  In 
the  latter. 

II.  EXPERIMENTAL 

To  determine  the  temperature  and  field  dependence  of  the  dielectric 
stiffness  above  the  Curie  point,  we  use  an  arrangement  similar  to  the 
one  used  by  Drougard  and  Young5,  as  shown  In  Fig.  1.  It  Is  basically 
a  capacitance  measurement  with  measuring  capacitance  much  larger  than 
the  sample  capacitance.  The  crystal  under  study  is  subject  at  the  same 
time  to  a  slowly  varying  DC  biasing  field  and  to  an  audio-frequency 
(8kHz)  sine  wave  of  constant  low  amplitude  (about  1/300  of  the  coercive 
field) . 

The  samples  were  single-crystal  plates  with  thickness  between  100 

2 

and  250  and  area  between  0.1  and  0.3  cm  .  Both  surfaces  perpendicular 

to  the  polar  axis  were  completely  covered  with  evaporated  metal  electrodes 
2 

to  avoid  clamping  .  Inside  the  oven,  the  samples  were  placed  In  a  closed 
copper  cylinder  filled  with  helium  and  freelv  suspended  by  their  silver 
lead  wires.  The  temperature  of  the  oven  was  controlled  bv  a  Leeds  &  Northrup 
Electromax  III  controller  which  provides  temperature  stabilization  to  within 
0.01  °C  for  a  period  of  20  minutes  or  more. 

In  the  paraelectrtc  phase,  the  polarization  of  the  sample  is  Induced 
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by  the  biasing  field  applied  to  the  electrodes.  At  each  temperature 

-1  2 

chosen,  the  X-Y  recorder  plots  the  curve  of  X  as  a  function  of  E  . 
From  Eqs.  (2)  and  (3)  we  get  an  expression  for  X  '  as  a  function  of 
E2. 


<~l  =  A  +  3B(E/A)2  +  (-6B2/A  +  5C)  (E/A)4  +  (21-~  -  26  ^  +  7D)  (|)6 

A 


By  least-squares  fitting  we  can  determine  the  coefficients  A,B  and  C. 

It  must  be  remembered,  however,  that  the  susceptibility  given  by 

Eqs.  (3)  and  (4)  is  the  Isothermal  one,  while  the  measurement  conditions 

2  6 

are  adiabatic.  It  is  easily  shown  ’  that  the  adiabatic  and  isothermal 
susceptibilities,  -v  and  Y»  respective • y,  are  related  by 


x^1  =  x^1  +  (t/b  cp)  (S2/*r)2  , 


(5) 


where  P  Is  the  density  and  Cp  Is  the  specific  heat  at  constant  pressure. 
With  the  help  of  Eq .  (2),  we  find 

XT1  =  A  +  3(B  +  TA  2/3p  C  )P2  +  5(C  +  2TA'b’/5p  C  )P4  +  ....  (6) 

3  P  p 


where  a  prime  Is  used  to  denote  differentiation  with  respect  to 
temperature.  The  usual  assumption,  that  only  A  depends  on  temperature, 
then  leads,  by  comparison  with  Eq .  (3),  to  a  correction  for  B, 


AB  =  TA  /3P  C  , 
o  p 


(7) 


(4) 
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and  no  corrections  for  the  higher-order  coefficients.  The  findings 
of  Ref.  2  and  the  predictions  of  Ref.  3  would  lead  to  corrections  in 
the  higher-order  coefficients  as  well,  e.g. 


AC  =  2TA  B  /5p  C  , 
oo  p 


(8) 


and  so  forth.  It  will  be  shown  below  that  the  corrections  for  C 
(and  presumably  the  higher-order  coefficients  also)  are  negligible  in 
the  materials  studied  here. 

III.  RESULTS 

A.  Barium  Titanate 

Equation  (4)  shows  that  only  near  the  Curie-Weiss  temperature 
T  ,  where  the  coefficient  A  approaches  zero,  can  we  expect  to  see 
large  effects  due  to  the  nonlinear  terms.  In  BaTiO^  a  first-order 

transition  occurs  at  a  temperature  significantly  above  Tq,  and  the 

-1  2 
measured  X  is  always  a  linear  function  of  E  ,  as  shown  in  Fig,  2. 

The  coefficient  A  has  the  expected  linear  behavior  (see  Fig.  3) 


with  the  values  of  A  and  T  listed  in  Table  I.  As  Indicated  in  the 

o  o 

table,  our  A  agrees  well  with  the  results  of  others,  but  T  is  10-15°C 
o  ’  o 

lower  than  generally  reported^.  Using  published  values  of  p  and 
in  Eq.  (7),  we  find  that  the  adiabatic  correction  is  approximately 


AB  *  1.8  x  10  cgs  throughout  the  temperature  range  studied,  Apply~ 


ing  this  correction  we  obtain  for  B  the  results  shown  in  Fig,  4,  The 

3  2 

linear  behavior  predicted  by  the  GMFT  and  reported  by  Drougard  et  al 


is  evident,  and  least-squares  fitting  gives  Bq  and  T^  as  listed  in 

2 

Table  I,  The  values  obtained  by  Drougard  et  al  are  given  for  comparison, 
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T 


B.  Triglycine  Sulfate 

Since  the  phase  transition  in  TGS  is  second  order,  measurements 

can  be  made  in  the  paraelectric  phase  quite  close  to  T  .  As  a  result, 

significant  nonlinearities  can  be  seen  when  x  ^  Is  plotted  as  a  function 
2 

of  E  ,  and  this  is  shown  in  Fig.  5.  For  temperatures  more  than  2  K  above 

Tq  the  relationship  is  essentially  linear,  but  closer  to  Tq  substantial 

4  6 

contributions  from  terms  of  order  E  and  E  are  required  to  fit  the 

experimental  data.  (The  solid  curves  in  Fig.  5  represent  least-squares 

fits  to  the  measured  values  indicated  by  crosses.)  The  fitting  gives 

A  immediately,  and  the  temperature  dependence  of  this  coefficient  is 

shown  in  Fig.  6.  The  linear  behavior  required  by  the  Curie-Weiss  law 

is  found,  and  the  values  of  Aq  and  T  ,  listed  in  Table  II,  are  in  good 

7  8 

agreement  with  the  results  of  others  ’  .  Again  making  use  of • published 
results  for  the  density  and  specific  heat,^  we  find  that  the  adiabatic 
corrections  for  B  are  approximately  10%  and  those  for  C  are  negligible. 
After  making  the  necessary  corrections  we  obtain  the  results  for  these 
coefficients  shown  in  Figs.  7  and  8.  In  each  case  the  temperature 
dependence  is  linear,  as  predicted  by  the  GMFT,  to  within  the  accuracy 
of  the  experiment.  Least-squares  fits  to  the  forms  B  =  Bq(T  -  T^)  and 
C  =  Cq(T  -  T2)  are  shown  'n  Figs.  7  and  8  and  the  fitting  parameters  list¬ 
ed  in  Table  II  which  also  gives  for  comparison  the  values  of  B  and  C 

g 

obtained  by  Triebwasser. 

IV.  DISCUSSION 

The  support  given  the  GMFT  by  the  results  of  the  previous  section 

is  substantial  but  not  complete.  In  the  base  of  BaTiO^,  our  findings 

differ  quantitatively  (but  not  qualitatively)  from  those  of  previous 
2 

studies  and  yield  somewhat  less  plausible  values  of  the  parameters  ap¬ 
pearing  in  the  GMFT. 
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3  4 

A  central  assumption  of  the  GMFT  *  is  the  replacement  of 

the  Lorentz  local  field  E^  *  E  -f  XP  by  a  nonlinear  local  field 

3 

correction  E,  -  E  +  X  P  +  X.P  +  ....  where  the  local  field  co- 
/  o  1 

efficients  are  given  by  X  ■  A  T  ,  X,  *  B  T, ,  and  so  forth.  Our 

o  o  o  1  o  1 

2 

experimental  results  and  those  of  Drougard  et  al  yield  similar 

values  of  XQ  and  X^,  which  are  listed  in  Table  III.  It  is  worth 

noting  that  these  products  of  experimentally  measured  quantities 

differ  less  between  the  two  sets  of  data  than  do  the  measured 

quantities  themselves.  In  addition,  the  values  of  X  calculated 

o 

9 

from  experimental  data  for  KNbO^  and  PbTiO^  differ  by  less  than 

25%  from  that  of  BaTiO^  despite  differences  of  as  much  as  a  factor 

of  two  in  T  and  A  .  Both  these  results  suggest  that  the  local 
o  o 

field  corrections  are  in  fact  characteristic  of  the  perovskite 
structure . 

The  GMFT  also  relates  the  coefficients  A  and  B  to  the  maxi- 

o  o 

mum  (saturation)  polarization, 

Np  =  (3Ao/5Bq)%  (9) 

and  the  density  of  independent  dipoles, 

N  =  A^/$kBo,  (10) 

where  k  is  Boltzmann's  constant.  Table  III  lists  the  values  of  these 

two  quantities  as  calculated  from  our  experimental  results  and  those 
2 

of  Drougard  et  al  .  (We  give  N  in  the  form  of  the  ratio  N  /N,  where 
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Nq  is  the  number  of  unit  cells  per  unit  volume,  since  this  is  thq 

number  of  unit  cells  required  to  form  one  of  the  elementary  dipoles  , 

assumed  by  the  GMFT.)  The  two  calculated  values  of  Np,  are  similar, 

2  2 

but  somewhat  smaller  than  the  26  \T./cm  obtained  experimentally 

Our  data  yield  a  significantly  larger  value  of  Nq/N  than  that  of 
2 

Drougard  et  al,  and  both  values  are  larger  than  the  estimate 

10  <  N  /N  <25  obtained  by  Lambert  and  Comes  from  measurements  of 
o 

the  diffuse  scattering  of  x-rays. ^  Estimates  based  on  the  observed 

entropy  change  at  the  transition  indicate,  however,  that  the  dipole 

length  is  highly  variable  from  sample  to  sample  and  may  be  as  great 

as  60  unit  cells. ^  Hence,  the  values  reported  in  Table  III,  while 

large,  are  by  no  means  unreasonable. 

In  the  case  of  TGS ,  we  find  that  each  of  the  coefficients 

measured  has  a  linear  temperature  dependence  (within  experimental 

accuracy)  throughout  the  temperature  range  investigated.  Unfortunately, 

however,  a  careful  examination  of  Eq .  (4)  and  our  experimental  data 

shows  that,  due  to  the  accidental  near  cancellation  of  several  terms 

in  the  calculation,  the  accuracy  with  which  the  coefficients  C  and  D 

4  6 

can  be  determined  from  the  coefficients  of  E  and  E  is  quite  poor. 

The  values  obtained  for  D  are  of  virtually  no  significance  and  will 
not  be  considered  here.  The  results  for  C  could  be  in  error  by  as 
much  as  a  factor  of  two.  Further  doubt  is  cast  on  the  values  obtained 
for  C  by  the  fact  that  using  the  coefficients  in  Table  II  one  calculates 
a  spontaneous  polarization  which  has  nearly  saturated  at  room  tempera¬ 
ture,  reaching  a  value  only  about  one  third  that  observed  experimentally. 
This  is  due  to  the  large  and  rapidly  increasing  value  of  C  implied  by 


the  linear  fit  shown  in  Fig.  8.  This  difficulty  would  disappear, 
however,  if  C  were  only  weakly  temperature  dependent.  Using  a 
constant  value  of  C  ■  1.9  x  10  ^  Cg8  we  obtain  excellent  agree¬ 
ment  between  the  calculated  spontaneous  polarization  and  that 
measured  by  Choe  et  al  ^ ,  as  shown  in  Fig.  9.  This  value 

g 

for  C,  although  larger  than  those  reported  by  Triebwasser  (see 

Table  II)  and  Gonzalo  ^ ,  is  close  to  that  found  by  Zerem  and 
12  -17 

Halperin  (1.4  x  10  cgs)  when  they  include  the  contribution 

of  the  coefficient  D  in  calculating  the  spontaneous  polarization. 

This  value  is  also  consistent  with  the  GMFT,  as  will  be  shown  below, 

and  is  well  within  the  range  of  our  experimental  results. 

Our  experimental  results  for  A  and  B  in  TGS,  when  substituted 

oo 

in  Eqs .  (9)  and  (10),  yield  strong  support  for  the  GMFT.  The  cal- 

2 

culated  saturation  polarization  is  Np  *  4.4  pc/cm  ,  in  excellent 

13  2 

agreement  with  the  measured  low- temperature  value  of  4.3  \fJ cm  . 
The  elementary  dipole  is  found  to  have  length  Nq/N  =»  1.1  unit  cells, 
in  full  accord  with  the  conventional  view  of  TGS  as  an  order-disorder 
ferroelectric . 

The  negative  value  found  for  CQ  (Table  II)  is  inconsistent  with 
the  GMFT,  however.  The  theory  predicts 

Co  *  99  Aq/175  (Np)\  (11) 

which  is  necessarily  positive.  Fortunately,  the  value  required 
above  to  obtain  good  agreement  with  the  observed  spontaneous  polari¬ 
zation  (C  ■  1.9  x  10  ^  cgs,  independent  of  temperature)  is  con- 
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sistent  with  both  Eq.  (11)  and  the  predicted  temperature  dependence, 

C  =  Cq(T  -  T2).  (12) 

If  the  value  quoted  for  C  Is  assumed  correct  at  T  =  Tq,  then  Eqs.  (11) 
and  (12)  can  be  solved  for  T2>  yielding  T2  =  46  K.  The  small  value 
of  T£»  when  used  in  Eq.  (12),  guarantees  that  C  is  weakly  temperature 
dependent,  varying  by  only  10%  over  the  temperature  range  in  which 
the  spontaneous  polarization  was  fitted  (Fig.  9). 

While  it  is  clearly  necessary  to  develop  more  accurate  methods 
of  determining  the  high-order  Devonshire  coefficients  in  order  to 
test  the  temperature  dependence  predicted  by  the  GMFT,  all  of  our 
results  are  consistent  with  that  theory,  and  the  most  reliable  of 
them  provide  strong  new  evidence  for  its  validity. 
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Table  I.  Devonshire  coefficients  of  BaTlO^  found  In  the  present  study 
and  those  reported  by  Drougard  et  al  (Ref.  2). 


Present  work  Reference  2 

A(10_5cgs)  7.3  (T  -  369.4K)  7.6  (T  -  378K) 

B(10-14  cgs )  2.3  (T  -  406K)  1.8  (T  -  448K) 
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Table  II. 


Devonshire  coefficients  of  TGS  found  In  the  present  study 
and  those  reported  by  Trlebwasser  (Ref.  8). 


Present  work 

Reference 

8 

A(10  3cgs) 

3.57  (T  -  321.8  K) 

3.92  (T  - 

323K) 

B(10  l°cgs) 

12.5  (T  -  277  K) 

8.0 

C(10"17  cgs) 

-8  (T  -  324  K) 

0.5 
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Table  III. 


Parameters  of  the  generalized  molecular  field  theory 


for  BaTiO^  (see  text)  as  calculated  from  the  results 
of  the  present  work  and  those  of  Drougard  et  al  (Ref.  2) 


Parameter 


Present  work 


Reference  2 


Fig.  1.  Apparatus  for  recording  the  dielectric  stiffness  as  a  function 
of  the  temperature  and  the  square  of  the  applied  electric  field. 
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(  1 /DIELECTRIC  CONSTANT 


Fig.  2.  The  dielectric  stiffness  x  of  BaTi03  as  a  function  of  the 
square  of  the  applied  electric  field  for  several  temperatures  above 
the  ferroelectric-paraelectric  transition. 
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TEMPERATURE  C 


Fig.  3.  The  first  Devonshire  coefficient  A  as  a  function  of  temperature 
in  BaTiO„ . 


(invisnoo  oiHioaiaia/ufr) 


Fig.  5.  The  dielectric  stiffness  X  of  TGS  as  a  function  of 
the  square  of  the  applied  electric  field  for  several  temperatures 
above  the  f erroelectric-paraelectric  transition. 
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TEMPERATURE  °C 
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The  first  Devonshire  coefficient  A  as  a  function  of  temperature 
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I.  INTRODUCTION 


It  ts  the  purpose  of  this  note  to  review  t!*e  known  limitations 
on  the  performance  of  ferroelectric  materials  as  pyroelectric  vidlcon 
targets  and  to  demonstrate,  both  theoretically  and  experimentally, 
that  a  relatively  new  class  of  materials,  the  Improper  ferroelec tries, 
can  avoid  these  limitations  and  may  eventually  yield  better  performance 
than  the  materials  now  in  use.  We  first  introduce  the  material  figure 
of  merit,  which  is  proportional  to  the  signal-to-noise  ratio  obtainable 
from  a  given  material  when  used  as  a  vidicon  target.  We  then  extend 
previous  results^  to  show  that  the  estimated  upper  limit  for  the  figure 
of  merit  originally  obtained  on  the  basis  of  a  particular  model  for 
ferroelectric  phase  transitions  is  in  fact  a  property  of  a  large  class 
of  generalized  molecular  field  models.  When  the  temperature  dependence 

of  the  figure  of  merit  is  examined,  using  both  the  generalized  molecular 

2  3 

field  theory  (GMFT)  and  the  Devonshire  formalism  ,  we  find  that  when 

either  of  these  theories  applies  the  performance  of  a  vidicon  target  can 
be  expected  to  deteriorate  as  its  Curie  temperature  is  approached,  despite 
a  large  increase  in  its  pyroelectric  coefficient.  Finally,  we  take  note 
of  a  class  of  materials,  the  improper  ferroelectrics,  to  which  these  con¬ 
siderations  do  not  apply  and  for  which  a  radically  different  temperature 

dependence  of  the  figure  of  merit  is  predicted.  Preliminary  experimental 
4 

results  confirm  this  theoretically  predicted  behavior  and  indicate  that 
a  well-known  improper  ferroelectric  (terbium  molybdate)  has  a  figure  of 
merit  near  its  transition  temperature  which  would  make  it  competitive  with 
the  best  of  the  materials  now  in  use  as  vidicon  targets. 


41 


II.  FIGURE  OF  MERIT 


It  Is  well  known^  that  the  signal-to-nolse  ratio  of  a  pyro¬ 
electric  vidicon  depends  on  the  properties  of  the  target  material 
through  the  ratio 


M  *  p/ec  (1) 

where  p  is  the  pyroelectric  coefficient,  e  is  the  dielectric  constant, 
and  c  is  the  volume  specific  heat.  In  the  following  treatment  of  the 
figure  of  merit  M,  we  will  not  consider  the  specific  heat  in  any  detail, 

since  its  variations,  both  with  temperature  and  from  one  material  to 

another,  are  small  compared  to  those  exhibited  by  the  pyroelectric  co¬ 
efficient  and  dielectric  constant.  It  is  also  convenient  to  replace 
the  dielectric  constant  by  the  susceptibility  X  *  e-l,  an  approximation 
which  is  quite  reasonable  in  view  of  the  large  values  of  e  typical  of 
ferroelectrics .  Using  these  approximations,  we  are  able  to  investigate 
the  behavior  of  the  figure  of  merit  M  by  studying  the  ratio 

m  *  p/x»  (2) 

which,  as  we  have  noted  previously^,  is  quite  easily  calculated  from 
most  models  for  ferroelectric  phase  transitions. 

III.  LOCAL  FIELD  MODELS 

We  have  reported  elsewhere*  that  the  existence  of  a  ferroelectric 

4  2 

crystal  with  a  value  of  M  significantly  greater  than  3.1  x  10  cm  / C 
is  most  Improbable.  (The  largest  value  reported  to  date^  is  that  of 


3  2 

trlglyeine  f luoberyllate ,  9.5  x  10  cm  /C .  )  This  estimated  upper 

limit  Is  based  on  consideration  of  a  local  (or  molecular)  field 

2  6  7 

model  which  has  been  described  In  detail  In  previous  reports.  ’  ’ 

We  can  now  show  that  a  large  class  of  local  field  models,  which  in¬ 
cludes  the  model  used  In  our  earlier  work,  yields  an  expression  for 
the  figure  of  merit  Identical  to  that  obtained  In  the  special  case 
studied  previously.  This  result  strengthens  our  belief  that  proper 
ferroelectrlcs  to  which  such  models  apply  offer  little  hope  for 
improved  pyroelectric  vldlcon  targets. 

In  a  local  field  model  of  a  ferroelectric  phase  transition  we 
represent  the  crystal  by  an  ensemble  of  N  electric  dipoles  per  unit 
volume,  where  each  dipole  has  a  moment  p,,  and  the  dipoles  interact 
via  a  local  field 


Ej ■  -  E  ♦  \(P),  (3) 

where  E  is  the  macroscopic  (applied)  field,  and  the  local  field 
correction  \(P)  is  an  arbitrary  function  of  the  polarization  P.  In 
order  to  apply  the  techniques  of  statistical  mechanics  to  such  a 
system,  it  is  necessary  to  specify  which  relative  orientations  of  the 
local  field  and  dipole  moment  are  permitted.  One  then  obtains  a 
relation  of  the  form 


P  =  Nn  f(x). 


(4) 


3 


where  x  »  pE./kT,  and  the  form  of  f(x)  is  determined  by  the  allowed 

V 

orientations  of  the  dipole  in  the  local  field.  Thus,  a  particular 
local  field  model  is  obtained  by  specifying  the  polarization  de¬ 
pendence  of  the  local  field,  i.e.  X(P),  and  the  permitted  dipole 

2  6  7 

orientations,  i.e.  f(x).  Our  previous  work  *  ’  has  dealt  with 

3 

local  fields  which  have  the  form  of  a  power  series  X(P)  =  XqP  +  X|P  + 
^P5  +.  . .  or  the  closely  related  Pade  approximant  A(P)  =  XoP/[l“(X^/Xo>P2 
and  with  classical  ("infinite  spin")  systems  in  which  the  dipoles 
can  have  any  orientation  relative  to  the  local  field,  so  that  f(x) 
is  the  Langevin  function  coth  x  -  l/x. 

For  an  arbitrary  local  field  model,  we  can  investigate  the  ratio 
m  «*  p/x  by  using  the  relation 

m  *  (dE/3T)p,  (5) 

which  follows  from  elementary  thermodynamics.^  (Note,  however,  that 
the  susceptibility  x  in  Eq.  (5)  is  the  isothermal  one  and  that  adiabatic 
corrections  may  then  be  needed  to  predict  the  performance  of  an  actual 
vidicon  target.^-)  Assuming  that  the  local  field  correction  X(P)  has 
no  explicit  temperature  dependence,  Eq*.  (3)-(5)  yield 

m  ■  kx/^.  (6) 

Equation  (6)  is  the  relation  which  was  obtained  previously  for  the 
CMFT  and  which  led  to  the  estimated  upper  limit  for  M  cited 

above.  We  now  see  that  this  result  is  more  general  than  originally 
supposed  and  applies  to  all  local  field  models. 


IV.  TEMPERATURE  DEPENDENCE  OF  THE  FIGURE  OF  MERIT 

The  general  result  (5)  can  also  be  used  to  study  the  tempera¬ 
ture  dependence  of  the  figure  of  merit.  If  we  wish  to  use  a  local 

field  model,  then  Eq .  (6)  applies  as  well.  In  each  of  the  local 

2  (y  "J  .  1 

field  models  studied  to  date,  *  *  we  have  found  that  x  =  f  (P/Np) 

is  a  rnonotonically  increasing  function  of  p,  which  vanishes  at 

P  =  0.  (We  expect  that  x  will  have  these  properties  in  all  physically 

reasonable  local  field  models,  and  work  is  under  way  to  explore 

these  and  other  general  characteristics  of  local  field  models.) 

Thus,  in  every  local  field  model  studied  so  far  the  figure  of  merit 

(or  more  precisely,  the  ratio  m)  is  found  to  decrease  as  the  Curie 

temperature  is  approached  and  the  spontaneous  polarization  decreases. 

A  similar  result  is  obtained  from  Devonshire 's  phenomenological 

3 

theory  of  ferroelectr icity .  If  the  free  energy  is  expanded  as  a 
power  series  in  the  polarization, 

ry  /  «i  r 

G  -  ’s  AP  +  i;  BP  +  £  CP  +. . .,  (7) 

then  the  electric  field,  E  =  (dG/SP)^,  is 

3  5 

E  ■  AP  +  BP  +  CP  -f. ..,  (8) 

and  Eq.  (5)  yields 

m  =•  A'P  +  B'P3  +  C'P5  +...,  (9) 
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+.K 


where  a  prime  indicates  differentiation  with  respect  to  temperature. 
It  is  usual,  in  applying  the  Devonshire  formalism,  to  assume  that 
A  has  the  form  Ao(T-To)  and  that  the  other  coefficients  are  tempera¬ 
ture  independent.  This  leads  to  the  simple  relation 


m  =  AqP.  (10) 

(The  corresponding  result  for  zero  applied  field,  m  =  AoFg,  where 
Pg  is  the  spontaneous  polarization,  has  been  obtained  elsewhere 
under  more  restrictive  assumptions.)  While  Eq.  (10)  clearly  pre¬ 
dicts  a  monotonic  decrease  in  m  at  zero  field  as  the  Curie  tempera¬ 
ture  is  approached,  the  question  of  the  temperature  dependence  of 

9 

the  higher-order  Devonshire  coefficients  is  by  no  means  resolved. 

The  available  experimental  evidence  indicates  that  these  coefficients 
are  in  fact  temperature  dependent,  but  that  the  observed  dependences 
are  at  least  consistent  with  the  prediction  of  the  GMFT  that  all  the 
coefficients  A1,  B*,  C1,...  are  positive,  so  that  Eq.  (9)  also  leads 
to  3  monotonically  decreasing  m. 

4 

Finally,  a  direct  measurement  of  the  temperature  dependence 
of  the  figure  of  merit  M  in  triglycine  sulfate  (TGS)  lends  strong 
support  to  our  theoretical  expectations.  Using  dynamic  response 
measurements  of  the  type  originated  by  Chynoweth  ^  and  further 
developed  by  Simhony  and  Shaulov,^  one  can  obtain  a  pyroelectric 
signal  (voltage)  which  is  directly  proportional  to  M.  This  technique, 
which  will  be  described  in  detail  elsewhere,  eliminates  the  need 
to  calculate  products  or  ratios  of  independent  measurements,  a  process 
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which  can  lead  to  serious  errors,  especially  in  temperature  ranges 
where  the  properties  involved  are  strongly  temperature  dependent. 

The  result  for  TGS  is  given  in  Fig.  1  and  clearly  shows  the  expected 
behavior.  The  small  rise  in  M  just  above  room  temperature  results 
from  a  reduction  in  the  pyroelectric  voltage  at  lower  temperatures 
due  to  loading  by  the  input  capacitance  of  the  measuring  amplifier 
and  is  completely  eliminated  when  one  corrects  the  data  for  this 
effect.  Since  the  coefficients  B*,  C',...  in  Eq.  (9)  are  not  known, 
measurements  of  this  type  can  yield  valuable  information  of  a 
fundamental  nature  as  well  as  predictions  of  vidicon  target  per¬ 
formance  . 

V.  IMPROPER  FERROELECTRICS  IN  VIDICON  APPLICATIONS 

The  considerations  of  the  previous  sections  indicate  that  among 
ferroelectrics  to  which  the  Devonshire  or  local  field  theories  apply, 
one  is  unlikely  to  obtain  a  material  with  a  significantly  larger 
room  temperature  figure  of  merit  than  those  now  in  use.  In  addition, 
attempts  to  improve  performance  by  operating  the  target  close  to  its 
Curie  temperature,  where  the  pyroelectric  coefficient  is  large,  must 
fail  due  to  the  temperature  dependence  of  M  predicted  by  Eqs.  (6)  and 
(9)  and  reflected  in  the  data  of  Fig.  1.  The  decrease  in  the  figure 
of  merit  as  the  Curie  temperature  is  approached  can  be  regarded  as  a 
consequence  of  the  fact  that  while  both  the  dielectric  constant  and 
pyroelectric  coefficient  Increase  dramatically  near  the  ferroelectr ic- 
to-pyroelectr ic  transition,  a  free  energy  of  the  form  (7)  will  always 
require  that  the  dielectric  constant  dominate  the  ratio.  However,  in 
recent  years  there  has  been  considerable  interest  in  an  unusual  class 


of  ferroelec Cries  which  cannot  be  described  by  the  free  energy  (7) 

and  which  typically  exhibit  only  a  small  increase  in  dielectric 

constant  near  the  phase  transition.  These  improper  or  extrinsic 

ferroelectrics  have  been  the  object  of  much  theoretical  and  ex- 
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perimental  investigation,  and  it  has  been  suggested  that, 
without  a  large  increase  in  dielectric  constant  to  compete  with 
the  increase  in  pyroelectric  coefficient,  these  materials  might 
exhibit  substantial  increases  in  figure  of  merit  as  the  transi¬ 
tion  temperature  is  approached.  In  the  following,  we  will 
demonstrate  both  theoretically  and  experimentally,  that  this 
suggestion  is  correct  and  could  lead  to  major  improvements  in 
pyroelectric  vid icon  performance. 

From  a  theoretical  point  of  view,  the  distinguishing  feature 
of  an  improper  ferroelectric  is  that  the  order  parameter  which 

describes  the  spontaneous  symmetry  change  at  the  paraelectr ic- 

14 

ferroelectric  phase  transition  is  neither  the  spontaneous 

polarization  nor  any  other  quantity  having  the  same  transforma- 
12 

tion  properties.  Instead,  a  polarization  appears  at  the  transi¬ 
tion  temperature  as  a  result  of  coup! ing  between  the  order  para¬ 
meter  (or  soft  mode)  and  the  polarization  (or  polar  optical  mode). 

This  coupling  can  be  described  in  macroscopic  terms  by  using  a 
free  energy  which  has  been  expanded  as  a  power  series  in  both  the 
order  parameter  11  and  the  polarization  P  and  which  contains  coupling 
terras  proportional  to  products  of  11  and  P.  The  simplest  such  free 
energy  has  been  treated  by  Dvorak  ^  and  has  the  form  (at  zero  stress) 


(11) 


G  =  h  A(T-To)  if  +  \  blf  +  i  c  7]6  +  ... 

+  3j  £  p2  +  fT)V 

The  electric  field  E  =  (dG/dP)^  is  given  by 

E  =  x"1  P  +  £lf,  (12) 

o 


and  applying  Eq .  (5)  yields 


m  -  nfTfl(dT)/9T)t 


(13) 
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The  derivative  in  Eq .  (13)  can  be  evaluated  by  noting  that  the 

)T 


spontaneous  value  of  the  order  parameter  7]  is  given  by  (dG/dll)T  =■  0 


or 


a (T-T  )  +  blf  +  elf  +  n'TT'2  P  =  0  . 
o  s  s  s 


(14) 


If  we  differentiate  Eq.  (14)  with  respect  to  T  and  substitute  in 
(13),  assuming  E  *  0,  we  obtain 

m  -  -nafT£‘2/[2(b+4cTl^)  -  n(n-2)  f2Xu^(n’2)]-  (15) 

Although  the  temperature  dependence  of  (15)  is  not  obvious,  3 
particularly  simple  case  is  provided  by  the  rare-earth  molybdates 
of  the  form  R2^Mo°4^3’  where  R  =  Gd>  Tb»  Sm>  9tc*  In  tbis  case 

u9 


Che  appropriate  free  energy  is  obtained  ^  for  n  =  2,  and  Eq .  (15) 
becomes 

m  =  -  af/(b+2c  ip.  (16) 

It  is  readily  shown  ^  that  the  stability  of  the  ferroelectric 
phase  requires  both  c  and  b  f-  2c  to  be  positive.  Hence,  m  has 

a  maximum  at  the  minimum  of  T]  ,  and  this  presumably  occurs  at  the 
Curie  temperature. 

Thus,  in  sharp  contrast  to  the  case  of  proper  ferroelectrics , 

we  expect  the  figure  of  merit  of  a  rare-earth  molybdate  to  increase 

as  the  transition  temperature  is  approached.  Figure  2  shows  M  as 
4 

measured  in  terbium  molybdate  (TMO)  by  the  dynamic  response 

technique.  A  detailed  discussion  of  the  results  for  TMO  and  several 

other  molybdates  will  be  given  elsewhere.  In  the  present  context  it 

is  sufficient  to  note  that  the  figure  of  merit  of  TMO  rises  from  a 

2 

room-temper ature  value  of  about  350  cm  /C  to  a  peak  at  the  Curie 

2 

temperature  greater  than  6000  cm  /C.  This  dramatic  increase  gives 
it  a  figure  of  merit  greater  (in  a  limited  temperature  range)  than 
that  of  TGS  and  approaching  that  of  the  best  available  material, 
triglycine  f luoberyllate .  The  rare-earth  molybdates  were  chosen 
for  this  preliminary  work  because  they  are  readily  available  and 
widely  studied.  No  attempt  was  made  to  select  or  modify  the 
material  to  achieve  a  large  figure  of  me -it.  For  this  reason  we 
believe  that  the  potential  of  improper  ferroelectrics  for  improved 
pyroelectric  vidicon  target  materials  is  largely  untapped  and  that 
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ABSTRACT 


A  study  of  the  published  literature  on  ferroelectric  thermal 
to  electrical  energy  conversion  provides  the  basis  for  classifying 
proposed  ferroelectric  converters  into  three  types,  based  on  their 
distinctly  different  principles  of  operation.  These  types  are 
designated  thermodielectric ,  pyroelectric,  and  pyro-ferroelectric, 
and  a  figure  of  merit  is  derived  for  each,  permitting  an  assess¬ 
ment  of  the  relative  merit  of  various  ferroelectric  materials  when 
used  in  converters  of  each  type.  Results  from  various  theories  of 
ferroelectricity  and  from  the  experimental  literature  are  used  to 
evaluate  and  compare  these  figures  of  merit,  and  an  attempt  is  made 
to  identify  promising  directions  for  future  studies  aimed  at 
achieving  efficient  thermal  to  electrical  energy  conversion  using 
ferroelectrics . 
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I.  Introduction 


The  advantages  of  using  ferroelectric  materials  for  the  con¬ 
version  of  heat  energy  directly  into  electrical  energy  are  that 
they  are  radiation  resistant,  they  can  produce  high  voltage  out¬ 
puts  so  that  there  is  no  need  for  up-conversion,  they  can  produce 
either  AC  or  DC  output,  they  have  flat  spectral  responses,  and  they 
can  have  very  high  specific  outputs  (i.e.  power  output  per  unit 
weight).  The  disadvantages  of  ferroelectric  converters  are  that 
they  require  a  periodic  heat  source  and  the  absolute  conversion 
efficiencies  are  low.  The  ferroelectric  converters  considered 
to  date  can  be  classified  into  three  types,  thermodielectric, 
pyroelectric,  and  pyro-ferroelectric,  according  to  their  different 
operating  principles. 

II.  Thermodielectric  Converter 

Heat  conversion  can  be  accomplished  with  any  charged  capacitor 

in  which  the  permittivity  of  the  dielectric  varies  with  temperature. 

A  typical  circuit  of  a  thermodielectric  converter  is  shown  in 

Fig.  1.  The  dielectric  material  in  the  parallel-plate  capacitor  is 

normally,  although  not  necessarily,  a  ferroelectric  material.  The 

capacitor  is  initially  charged  by  the  storage  battery  through  the 

charging  diode.  At  the  Curie  temperature  Tc  the  capacitance  is  a 

maximum  because  of  the  large  dielectric  constant  e^.  The  charge  on 

the  capacitor  is  Q  »  CV,  where  V  is  the  storage  battery  voltage. 

If  the  capacitor  is  cooled  to  T  or  heated  to  T,  ,  the  dielectric 

a  b 

constant  drops  to  e  or  e  ,  and  the  capacitance  drops  by  a  factor 
a  d 

e^/e^  or  e^/e^.  The  excess  charge  which  appears  at  the  capacitor 
due  to  the  drop  in  capacitance  will  be  discharged  through  the 
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discharge  diode  and  load  the  storage  battery,  thus  completing  a 

thermal  to  electrical  energy  conversion  cycle. 

In  1959,  S.R.  Hoh  first  established  the  feasibility  of  thermo- 

1-3 

dielectric  energy  converters  using  ferroelectric  materials  .  In 

4 

1961,  Clingman  and  Moore  calculated  the  theoretical  conversion  . 
efficiency  for  polycrystalline  BaTlO^  and  found  that  it  is  in  the 
range  of  0.5%  to  1%.  Childress  ^  considered  the  sharper  phase  transi¬ 
tions  of  poled  single  crystals  of  BaTlO^,  included  the  field  induced 
phase  transition,  and  found  an  efficiency  of  0.5%.  In  1963,  S.R.  Hoh  ^ 
reviewed  the  subject  and,  instead  of  using  the  parameters  of  BaTjLO^, 
chose  a  set  of  idealized  material  parameters  and  found  an  efficiency 
of  5.6%.  Hoh  also  did  some  experimental  work  and  tried  to  use  thin 

films,  but  was  unable  to  prepare  suitable  ferroelectric  films.  In 
7  8 

1963,  L.  Eyraud  and  J.  Perrigot  in  France  also  studied  the  thermo- 

9 

dielectric  converter.  Recently,  in  1976,  Pulvari  and  Garcia  proposed 

the  use  of  a  multiple-layer  structure  to  achieve  higher  dielectric 

breakdown  and  possibly  a  higher  conversion  efficiency. 

4 

Analysis  of  the  efficiency  of  the  thermodielectric  converter 

? 

leads  to  consideration  of  the  ratio  CE“/c,  where  C  is  the  Curie  constant 

D 

E^  is  the  dielectric  strength  (breakdown  field),  and  c  is  the  volume 

specific  heat.  The  larger  this  figure  of  merit,  the  better  a  material 

will  perform  in  a  thermodielectric  converter.  Using  published  values 

4 

for  the  parameters  of  BaTiO^,  one  finds  that  the  figure  of  merit  is 

15  2  2  3 

3.1  x  10  V  cmK  /cal.  A  realistic  efficiency  for  a  BaTiO^  thermo- 

-3 

dielectric  converter  is  about  10  ;  a  practical  converter  would  require 

a  material  with  a  figure  of  merit  at  least  an  order  of  magnitude  greater 
than  that  of  BaTiO^. 
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III.  Pyroelectric  Converter 

A  pyroelectric  converter  (Fig.  2)  uses  the  pyroelectric  effect  to  gener¬ 
ate  electrical  power.  The  incoming  light  is  chopped  and  the  variation 
in  sample  temperature  is  small,  with  the  temperature  of  the  material 
remaining  below  the  transition  temperature.  This  technique  was  con¬ 
sidered  by  L.J.  Sivian  ^  as  early  as  1942,  and  in  1974  A.  van  der  Ziel 
showed  by  using  a  small-signal  model,  that  the  conversion  efficienty 
is  small  (0.025%)  for  most  pyroelectric  materials.  From  van  der  Zlel's 
derivation  we  see  that  the  influence  of  material  parameters  on  the  con- 
version  efficiency  is  proportional  to  p  /ce,  where  p  is  the  pyroelectric 
coefficient,  c  is  the  volume  specific  heat,  and  e  is  the  dielectric  constant. 
We  shall  use  this  quantity,  p  /ce,  as  the  application  figure  of  merit  for 
pyroelectric  energy  conversion  materials. 

IV.  Pyro-ferroelectric  Converter 

The  pyro-ferroelectric  converter  operates  with  the  same  principle 

as  the  pyroelectric  converter  except  that  the  variation  in  temperature 

is  large,  and  the  ferroelectric  material  is  heated  very  close  to  the 

transition  temperature  Tc ,  but  kept  slightly  below  to  avoid  depoling.  This 

achieves  maximum  utilization  of  the  spontaneous  polarization.  Because 

the  variation  in  temperature  is  large,  the  small-signal  theory  as  was 

used  by  van  der  Ziel  in  calculating  the  efficienty  of  pyroelectric 

12 

converter  is  not  valid  in  this  case.  Gonzalo  analyzed  the  stored 
and  the  convertible  energy  of  ferroelectrics  from  the  point  of  view 
of  their  basic  electrostatic,  thermal  and  transport  properties,  and  he 
concluded  that  ideal  efficiencies  of  the  order  of  10  -  15%  could  be 
attained  if  proper  materials  are  developed  to  be  used  in  multistage 
pyro-ferroelectric  converters. 
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According  to  Gonzalo,  the  conversion  efficiency  of  a  single 

stage  pyro-ferroelectric  converter  using  second  (or  close  to  second) 

order  proper  ferroelectrics ,  neglecting  strain  energy,  is  H  =  2rrp g/Cc 

(cgs  units)  where  C  is  the  Curie  constant,  c  is  the  average  volume 

specific  heat  over  AT,  and  Pg  is  the  spontaneous  polarization  at  the 

2 

temperature  Tc  -  AT,  with  the  requirement  that  dPg/dT  is  constant 

2 

within  AT-  Obviously,  we  can  use  the  quantity  Pg/Cc,  which  is  equal 
to  T1 /2tt,  as  the  application  figure  of  merit  for  proper  ferroelectrics 
used  as  pyro-ferroelectric  energy  conversion  materials.  For  the 
material  NaNO^,  with  AT  =  5  K  Gonzalo  found  a  conversion  efficiency  of 
0.27„  and  a  specific  output  of  4.4  x  10^  W/lb  (if  AT  =  75  K  then  one 
obtains  37,  and  1.5  x  10^  If/lb). 

« 

V.  New  Directions  and  Materials 

2 

The  figure  of  merit  p  /ce  which  determines  the  efficiency  of  a 
pyroelectric  converter  (Section  Tit)  might  be  expected  to  increase 
at  temperatures  near  a  ferroelectric-to-paruelectric  phase  transi¬ 
tion  where  the  pyroelectric  coefficient  p  increases  sharply.  Un¬ 
fortunately,  in  most  ferroelectrics  the  dielectric  constant  e  also 

13 

Increases  near  the  phase  transition,  and  it  can  he  shown  that 

14 

in  proper  ferroelectrics,  to  which  the  usual  Devonshire  theory 
applies,  the  ratio  p/(e-l)  is  proportional  to  the  spontaneous 
polarization  Pg  and  so  is  a  decreasing  function  of  temperature. 

In  recent  years,  another  class  of  ferroelectric  phase  transitions 
|  has  been  identified,  however,  in  which  the  polarization  is  not 

the  order  parameter. Such  improper  ferroelectrics  generally 
,  have  a  dielectric  constant  which  is  weakly  temper it-are  dependent. 
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often  exhibiting  only  a  small  discontinuity  at  the  Curie  point. 

As  we  have  noted  elsewhere,^3  this  behavior  can  lead  to  large 
values  of  p/e  near  th  j  transition,  and  recent  measurements  ^  on 
several  rare-earth  molybdates  have  confirmed  this.  Typical  results 
(for  gadolinium  molybdate,  Qd^CMoO^)^,  a  few  tenths  of  a  degree 
below  T  )  yield  p'/ec  -  l.l  x  10  "  / J  craK,  a  value  nearly 

2.5  times  greater  than  used  by  van  ler  2iel  in  hia  efficiency^ 
estimate.^  Obviously,  this  figure  of  merit  is  still  far  too 
small  to  produce  an  efficient  converter,  but  the  rare-earth  molybdates 
comprise  only  the  first  of  many  families  of  improper  ferroelectrics 
whose  unusual  properties  are  still  under  investigation,  and  there  is 
as  yet  no  reason  to  believe  that  very  much  greater  values  of  p  7ec 
will  not  he  found. 

A  second  remarkable  property  of  certain  improper  ferroelectrics 
forms  the  basis  of  a  novel  energy  converter  which  won 1 1  not,  in 
general,  be  feasible  using  proper  ferroelectrics.  We  have  recently 
observed  ^  that  despite  repeated  cycling  through  the  ferroelectric- 
paraelectric  phase  transition,  3n  essentially  single-domain  sample 
of  a  rare-earth  molybdate  will  retain  its  original  direction  and 
magnitude  of  spontaneous  polarization.  Levanyuk  and  Sannikov 
have  pointed  out  on  theoretical  grounds  that  for  an  improper  ferro¬ 
electric  such  transitions  from  a  parse lectric  to  a  single-domain 
ferroelectric  state  need  not  he  as  energetically  unfavorable  as  in 
a  proper  ferroelectric.  A  small  amount  of  elastic  energy,  perhaps 
due  to  clamping  or  internal  strains,  might  be  sufficient  to  stabilize 
the  single-domain  state.  With  the  possibility  of  depoling  eliminated 
In  thi.3  way,  we  may  consider  energy  conversion  cycles  in  which  the 
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temperature  passes  through  the  transition  temperature.  Only  the 
simplest  and  most  obvious  such  cycle  will  be  considered  here.  Re¬ 
finements  are  likely  to  lead  to  improved  efficiency. 


Let  us  suppose  that  a  slab  of  improper  ferroelectric  in  the 

form  of  a  parallel-plate  capacitor  with  area  A  and  thickness  d  is 

heated  through  a  first-order  phase  transition  by  an  incident 

radiation  flux  F,  and  then  cooled  back  through  the  transition  by 

interrupting  the  flux.  If  this  process  is  repeated  continuously, 

1 8 

it  is  readily  shown  that  when  steady  state  conditions  are 
achieved,  the  time  required  for  a  complete  heating  and  cooling 
cycle  is  r  =  4dA0/F,  where  AQ  is  the  net  heat  absorbed  (emitted) 
per  unit  volume  during  the  heating  (cooling)  part  of  the  cycle. 

A  short-circuit  current  of  average  magnitude  I  =  2PgcA/T  would 
flow  during  each  cycle,  where  Pgc  is  the  (discontinuous)  change  in 
Pg  at  Tc,  and  the  factor  of  two  arises  from  the  fact  that  equal 
but  opposite  currents  are  produced  during  the  heating  and  cooling 
periods.  If  a  load  resistor  R  is  placed  in  parallel  with  the 
sampl  !,  it  is  easily  shown  ^  that  the  power  dissipated  in  the 
load  is  a  maximum  for  R  =  t/2tC,  where  C  is  the  average  value  of 
the  sample  capacitance  during  the  cycle.  If  we  approximate  the 
current  in  the  load  as  a  square  wave  and  use  the  relation 
C  *  «A/4nd,  where  e  is  the  average  dielectric  constant,  then  we 
find  for  the  converter  efficiency  1]  =  I  R/FA  =  2Pgc/eAQ.  Assuming 
that  the  converter  is  operated  in  a  very  small  temperature  range 
in  the  vicinity  of  Tc ,  we  may  identify  AQ  with  the  latent  heat  of 
the  phase  transition  and  take  for  e  the  average  of  the  limiting 
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values  of  s  as  T£  is  approached  from  above  and  below. 

In  the  case  of  Gd2(MoO^)g,  the  material  whose  behavior  suggested 

this  conversion  cycle,  published  values  of  the  dielectric  and  thermal 

properties  ^  yield  an  efficiency  7)  =  0.4%.  Other  improper  ferro- 

electrics  offer  the  possibility  of  higher  efficiencies,  although 

their  tendency  to  enter  the  ferroelectric  phase  in  a  single-domain 

state  has  not  been  verified.  Iron-iodine  boracite  (Fe^B^O^I),  for 

example,  yields  an  efficiency  of  7]  ~  707.  when  the  data  of  Schmid  et 
19 

al  are  used.  These  results  indicate  that  as  investigation  con¬ 
tinues  into  the  behavior  of  the  boracites  and  other  families  of 
improper  ferroelectr ics ,  efforts  should  be  made  to  determine  their 
usefulness  in  converters  of  the  type  described  here. 

One  final  point  is  worth  noting.  The  conversion  cycle  proposed 
above  could  be  used  with  proper  as  well  as  improper  ferroelect"ics  if 
a  means  could  be  found  to  establish  a  preferred  domain  orientation 

as  the  material  enters  the  ferroelectric  phase.  Work  on  lanthanum 

20 

doped  lead  zirconate  titanate  ceramics  has  demonstrated  that 

space-charge  effects  in  aged  ceramics  can  establish  just  such  a 

preference,  and  this  possibility  is  now  under  investigation.  In  the 

case  of  a  proper  ferroelectric  with  a  first-order  transition,  the 

efficiency  obtained  above  reduces  to  71  =  8(Tc-Tq) /5tiTc  .  Using  the 

21 

values  tabulated  by  Jona  and  Shirane,  we  find  efficiencies  for 
BaTiO^  and  KNbO^  of  1%  and  4%,  respectively.  Although  these  effici¬ 
encies  must  be  reduced  to  account  for  the  fact  that  the  remanent 
polarization  of  a  ceramic  is  less  than  the  spontaneous  polarization 
of  the  corresponding  single  crystal,  they  are  large  enough  to  suggest 
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Anomalous  scattering  and  asymmetrical  line  shapes  in  Raman  spectra 
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The  Raman  spectrum  of  orthorhombic  KNbOj  is  found  to  contain  anomalous  scattering,  consisting  of  a  large 
background  and  broad  bands,  which  can  be  related  directly  to  the  dynamic  disorder  which  has  been  invoked 
to  explain  the  results  of  earlier  studies  of  diffuse  x-ray  scattering  and  inelastic  neutron  scattering.  Some 
ordinary  first-order  lines  exhibit  coupling  to  the  anomalous  part  of  the  spectrum,  which  is  revealed  by  the 
existence  of  resonant  interference  of  the  type  described  by  Fano.  AH  these  features  disappear  abruptly  at  the 
transition  to  the  low-temperature  rhombohedral  phase,  indicating  that  they  are  characteristic  of  the  peculiar 
dynamics  of  linear  chains  observed  in  the  orthorhombic,  tetragonal,  and  cubic  phases. 


INTRODUCTION  EXPERIM  ENTAL 

In  the  ferroelectric  perovskites  such  as  BaTiO,  The  Raman  scattering  was  excited  by  an  argon- 

and  KNbO,,  x-ray  diffuse-scattering1  and  inelas-  ion  laser  operating  at  either  4880  or  5145  A.  Mea- 

tic -neutron -scattering2-1  experiments  yield  evi-  surements  were  made  using  a  J-m  Spex  double- 

dence  of  a  very  large  anisotropy  in  the  dispersion  grating  spectrometer  located  at  the  Maybaum  In- 

of  some  of  the  vibrational  modes.  Light  scat-  stitute  of  Yeshiva  University  and  a  Coderg  T800 

tering','i  has  been  used  to  investigate  the  zone-  triple-grating  spectrometer  located  at  D^parte- 

center  phonons  in  order  to  check  the  soft-mode  ment  de  Recherches  Physiques,  University  Paris 

theories.'1  We  report  here  some  new  results  of  VI.  Both  systems  employed  photon-counting  de- 

Raman  scattering  experiments  performed  on  a  tection. 

single-domain  KNbO.,  sample,  both  in  its  room-  We  used  a  single-domain  sample  cut  with  the 

temperature  orthorhombic  phase  and  as  the  tern-  faces  perpendicular  to  the  orthorhombic  axes, 

perature  is  lowered  through  the  orthorhombic-to-  Standard  right-angle  and  backscattering  geom- 

rhombohedral  phase  transition.  etries  were  used  to  obtain  various  linear  momen- 

The  assignment  of  the  ordinary  first-order  lines  turn  transfers  from  the  photons  to  the  crystals, 

is  quite  easily  obtained,  leaving  outside  this  clas-  Table  I  summarizes  the  modes  which  can  be  ob- 
sification  a  large  continuous  scattering  which  we  served  in  each  scattering  geometry.  The  ortho- 

refer  to  as  “anomalous.”  The  latter,  previously  rhombic  axes  are  related  to  the  pseudocubic  ones 

described  in  Ref.  5  for  restricted  scattering  ge-  as  described  in  Fig.  1.  At  zero  wave  vector,  the 

ometries,  has  now  been  found  to  occur  more  gen-  irreducible  representations  of  the  12  optic  modes 

erally  whenever  the  incident  and  scattered  beams  in  the  C2„  point  group  of  the  orthorhombic  phase 

propagate  perpendicular  to  the  Z  axis  and  are  are  as  follows:  4A,(Y),  4B,(F),  3B2(Z),  1A,.  All 

simultaneously  polarized  parallel  to  ^.  As  ex-  An  B,,  and  B.4  modes  are  infrared  active  (the  let- 

pected,  this  anomalous  spectrum  disappears  ters  in  parentheses  indicate  the  direction  of  the 

abruptly  at  the  orthorhombic-to-rhombohedral  dipole  moment  of  each  mode),  and  for  these  the 

transition  temperature,  as  does  the  normal  low-  longitudinal  or  transverse  character  must  be 

frequency  8/ TO)  line.  In  the  spectra  where  the  considered.  Two  features  of  this  table  should  be 

anomalous  scattering  occurs,  some  first-order  noted.  First,  B,  and  B2  modes  are  not  detected 

lines  display  an  asymmetrical  shape  character-  in  a  longitudinal  configuration.  Second,  when  a 

istic  of  a  “Fano  interference.”7  polar  mode  is  neither  transverse  nor  longitudinal 
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TABLE  I.  Summary  of  the  scattering  geometries  employed  and  the  mode  symmetries  observed.  Scattering  genmolrie.s 
are  indicated  by  the  conventional  notation  where  kj.c,  (k,,?,)  are  the  wave  vector  and  polarization  '.-.  tor  ot 

the  incident  (scattered)  photon. 


Incident 

wave  vector 

Scattered 

wave  vector 

Momentum 

transfer 

Polarization  oi  incident  and  scattered  bo: 

and  mode  symmetry 

XX  YY  ZZ  YX.XY  ZX.XZ 

ims  U,?s' 

ri  ,y.) 

♦ 

k> 

k4-ki 

A, 

A, 

A, 

B, 

A  2 

Right-angle 

X 

Y 

±z 

.  .  . 

Mixed 

Mixed 

TO 

Indiff. 

scattering 

Y 

Z 

XX 

TO 

Mixed 

Mixed 

Indiff. 

X 

Z 

x  Y 

.  .  . 

Mixed 

.  .  . 

TO 

Mixed 

Indiff. 

X 

X 

II* 

... 

LO 

LO 

Indiff 

Backscattering 

Y 

? 

II Y 

TO 

TO 

Z 

2 

II  z 

TO 

TO 

. . . 

TO 

in  a  given  scattering  geometry  (labeled  “mixed”), 
it  loses  its  zero-wave-vector  irreducible-repre¬ 
sentation  species  through  mixing  with  modes  of 
different  symmetry. 

ROOM  TEMPERATURE  SPECTRA 

All  the  spectra  listed  in  Table  1  have  been  re¬ 
corded.  Only  those  of  special  interest  are  re¬ 
produced  in  Figs.  2  and  3.  In  Fig.  2,  the  peaks 
which  arise  from  contamination  by  other  species 
are  drawn  with  broken  Vines. 

From  Table  I,  and  taking  into  account  the  pos¬ 
sible  mixing  of  species,  we  can  assign  the  ob¬ 
served  frequencies  to  the  different  modes.  These 
assignments  are  summarized  in  Table  II.  Within 
the  limit  of  experimental  accuracy,  we  agree  with 
the  results  published  independently  by  Winter 
el  al.  ,5  at  least  for  those  lines  identified  unambig¬ 
uously  in  Ref.  5.  The  following  details  should  be 
noted.  The  low-frequency  TO)  line  (maximum 
at  40  cm'1)  is  broad  and  asymmetrical  [Fig.  2(a)], 
and  a  thin  line  at  195  cm'1  which  sits  on  its  tail  is 


FIG.  1.  Orthorhombic  axes  for  KNbOj  and  corre¬ 
sponding  pseudocublc  axes  (In  brackets).  Also  shown  are 
the  symmetry  elements  of  the  C2v  point  group.  The 
spontaneous  polarization  is  In  the  X  direction. 


slightly  distorted  and  seems  to  show  a  dip  on  the 
low-frequency  side.  In  the  Bz  mixed  spectrum, 
the  lines  at  170  and  205  cm'1  display  an  interfer¬ 
ence  dip  between  them  [Fig.  2(b)]. 

The  anomalous  spectrum  is  seen  in  the  following 
geometries:  Y(ZZ)Y,  X(ZZ)X,  and  X(ZZ)Y. 

These  are  A,(TO),  A,(LO),  and  Admixed)  spectra, 
respectively,  and  are  shown  in  Figs.  3(e),  3(a), 
and  3(c).  In  each  case,  the  linear  momentum 
transfer  is  perpendicular  to  the  Z  axis,  and  the 
light  beams  are  polarized  parallel  to  Z.  The 
anomalous  scattering,  identical  in  the  three  spec¬ 
tra,  consists  of  a  continuous  background,  de¬ 
creasing  in  intensity  from  low  to  high  frequency, 
and  two  broad  bumps  centered  near  130  and  430 
cm'1. 

The  ordinary  first-order  lines  are  superimposed 
on  the  anomalous  spectrum,  but  some  of  them 
show  a  coupling  leading  to  an  asymmetrical  line- 
shape  and  interference  dip.  This  is  clearly  dem¬ 
onstrated  by  comparison  of  the  pairs  of  spectra 
in  Fig.  3: 

(i)  Figures  3(a)  and  3(b):  The  normal  X(YY)X 


600  600  40O  300  100 

RAMAN  SHIFT  (cm'1) 


FIG.  2.  B 2  symmetry  spectra:  (a)  flj(TO)  in  back- 
srattering  geometry,  (b)  B2  (mixed)  in  right-angle  geom¬ 
etry.  The  dashed  lines  indicate  contamination  from 
other  symmetry  species  due  cither  to  a  small  misort- 
entation  of  the  polarizers  or  to  mixing  of  polar  modes  ns 
described  in  the  text. 


68 


5070 


QUIT  III,  B  K  L  L  ,  KRAUZMAN,  AND  RATfAH 


14 


FIG.  .7.  Anomalous  spectra  and  the  associated  normal 
spectra.  In  the  K(XX)F  and  Y (7.7.)'?  spectra,  the  broad 
line  in  the  vicinity  of  50  cm*  1  is  most  probably  a  con¬ 
tamination  from  the  Y(X7.)Y  B;(TO)  spectrum  due  to  the 
divergence  of  the  scattered  beam. 


spectrum  together  with  the  anomalous  X(ZZ)X 
one  are  both  A,(LO)  spectra.  Note  the  lines  at 
195  and  435  cm'1. 

(ii)  Figures  3(d)  and  3(e):  The  normal  Y(XX)Y 
spectrum  together  with  the  anomalous  X( ZZ)  Y 
one  are  both  A,(TO)  spectra.  Note  the  line  at 
193  cm'1. 

(iii)  Figure  3(c):  The  anomalous  X(ZZ)Y  spec¬ 
trum  has  A1  (mixed)  symmetry.  Note  the  mode  at 


TABLE  11.  Frequencies  and  symmetry  assignments  of 
first-order  Kaman  lines  of  KNbOj.  The  figure  in  paren¬ 
theses  is  the  full  width  of  the  lino  at  half-maximum. 


A,  (TO) 

A  ,(1,0) 

B  ((TO) 

B,(TO) 

A , 

40(27x2) 

197(2) 

194.5(3) 

192(2) 

249(25) 

270 

1.90.5* 

281.5(37) 

295(5) 

282(6) 

297(5) 

434.5*(10) 

534(20) 

513(20) 

934(27) 

“Lines  showing  an  asymmetrical  shape  in  at  least  one 
scattering  geometry. 


FIG.  4.  Temperature  dependence  of  part  of  the  B, 

+  B2  spectrum:  (a)  in  the  orthorhombic  phase  just  above 
the  phase  transition  (temperature  decreasing),  (b)  in 
the  rhombohedral  phase  just  below  the  transition  (tem¬ 
perature  decreasing),  (e)  in  the  rhombohedral  phase  in 
the  hysteresis  interval  (temperature  increasing). 

420  cm'1  which  appears  almost  wholly  negative  in 
the  anomalous  part  of  the  spectrum. 

TEMPERATURE  DEPENDENCE  (REF.  81(23  to  -70 °C) 
As  the  temperature  is  lowered,  there  is  essen¬ 
tially  no  change  in  the  low-frequency  B2(TO)  line 
and  the  anomalous  spectrum  until  the  transition 
temperature  is  reached.  The  frequency  and  width 
of  the  B2  line  both  decrease  slightly,  but  its  in¬ 
tensity  |  Fig.  4(a)  1,  as  well  as  the  shape  and  in¬ 
tensity  of  the  anomalous  spectrum,  remains  un¬ 
changed.  Below  the  phase  transition,  which  takes 
place  near  _59°C,  the  spectrum  becomes  very 
different  [Figs.  4(b)  and  4(c) [  as  the  broad  B.z  line 
and  the  anomalous  scattering  disappear  com¬ 
pletely.  Since  the  sample  becomes  broken  and 
twinned,  all  the  spectra  are  superimposed,  so 
that  it  is  not  possible  to  miss  any  lines.  Upon 
reheating  the  sample,  the  anomalous  spectra  of 
the  orthorhombic  phase  appear  again  at  -30  "C, 
exhibiting  the  large  thermal  hysteresis  character¬ 
istic  of  a  first-order  transition.  These  experi¬ 
ments  thus  demonstrate  that  the  anomalous  spec¬ 
trum  is  related  to  an  intrinsic  property  of  the 
orthorhombic  phase,  as  is  the  low-frequency 
B2(TO)  line. 


DISCUSSION 

We  note  first  that  all  the  anomalous  features 
which  appear  in  Raman,  neutron,  and  x-ray  scat¬ 
tering  experiments  occur  only  when  the  linear  mo¬ 
mentum  transfer  lies  perpendicular  to  the  Z  axis, 
and  that  these  features  disappear  in  the  rhombo¬ 
hedral  phase.  Hence  the  dynamical  properties  of 
the  orthorhombic  phase  have  a  pronounced  one¬ 
dimensional  character.  These  results  can  be 
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KIG.  5.  (a)  B2  spectrum  for  a  momentum  transfer  ro¬ 
tated  bv  15°  from  a  TO  configuration  toward  the  LO  con¬ 
figuration.  (b)  Anomalous  spectrum  for  the  same  mo¬ 
mentum  transfer. 


understood1  as  arising  from  displacements  of  the 
niobium  atoms  in  the  Z  direction  (away  from  the 
C,  symmetry  axis)  which  are  correlated  for  some 
finite  distance  along  Z.  Averaged  over  distances 
large  compared  to  this  correlation  length,  how¬ 
ever,  the  structure  retains  orthorhombic  sym¬ 
metry.  Such  disorder  need  not  be  static,  since 
its  effects  will  be  seen  in  scattering  experiments 
provided  the  lifetime  of  the  local  departure  from 
the  average  symmetry  is  greater  than  the  scat¬ 
tering  time.  From  this  point  of  view,  the  anom¬ 
alous  features  in  the  Raman  spectra  can  be  in¬ 
terpreted  as  disorder-induced  first-order  scat¬ 
tering,  i.e. ,  one-phonon  scattering  by  modes 
which  would  not  be  Raman  active  in  a  perfect 
crystal,  but  which  become  so  when  both  the  trans¬ 
lational  and  point-group  symmetries  are  broken 
by  the  disorder.  Since  all  phonons  (not  only  those 
at  the  zone  center)  can  become  Raman  active  in 
this  way,  one  anticipates  a  scattering  spectrum 
proportional,  in  first  approximation,  to  the  one- 
phonon  density  of  states.  This  interpretation  is 
further  strengthened  by  the  fact  that  neutron  scat¬ 
tering  results3  suggest  the  existence  of  a  peak  in 
the  one-phonon  density  of  states  near  150  cm'1, 
close  to  one  of  the  peaks  (130  cm'1)  in  the  anom¬ 
alous  Raman  spectra. 

We  have  studied  the  modification  of  the  spectra 
which  occurs  as  the  linear  momentum  transfer  is 
rotated  from  the  V  direction  to  the  Z  direction. 

As  the  momentum  transfer  is  changed,  the  inten¬ 
sity  of  the  anomalous  scattering  decreases  rapidly, 
while  the  frequency  of  the  B2  line  shifts  very  rap¬ 
idly  toward  higher  values  (Figs.  5(a)  and  5(b) |. 

We  conclude  from  this  that  the  B2  line  is  produced 
by  an  underdamped  zone -center  phonon  with  a 
large  LO-TO  frequency  difference  and  that  the 
anomalous  spectrum  in  ZZ  polarization  arises 
from  a  continuum  of  states.  Our  identification  of 


the  line  at  40  cm'1  as  the  zone-center  B2(TO)  pho¬ 
non  is  consistent  with  the  results  of  neutron  scat- 
tering'experiments.3  Although  Currat  cl  at.1  ex¬ 
trapolate  the  frequency  of  the  optical  branch  to 
about  25  cm'1  at  the  zone  center,  the  experi¬ 
mental  resolution  is  such  that  spectra  near  the 
zone  center  and  at  low  frequencies  are  always 
heavily  contaminated  by  the  Bragg  peak  and  the 
acoustic  dispersion.  Moreover,  the  experimental 
points  closest  to  the  zone  center  could  only  be  ob¬ 
tained  from  constant-energy  scans  which  do  not 
provide  an  accurate  measure  of  mode  frequency. 

The  asymmetrical  line  shapes  observed  in  this 
work  are  of  two  types.  One  is  produced  by  the 
small  interference  between  the  two  lowest-lying 
lines  in  the  B2(TO)  spectrum  [Fig.  2(a)  |.  This  is 
merely  a  coupling  between  modes  of  identical  sym¬ 
metry,  which  appears  because  of  the  large  width 
of  the  low-frequency  line.  The  coupling  becomes 
much  stronger  when  the  wave  vector  is  rotated 
45°  from  the  Y  axis  toward  the  Z  axis,  since  the 
rapid  increase  in  the  frequency  of  the  lower  mode 
brings  the  two  modes  quite  close  | Fig.  2(b)  |. 

The  second  type  of  asymmetry  can  be  seen  in 
the  shape  of  several  narrow  lines  observed  in  the 
anomalous  ZZ  spectra.  This  can  be  analyzed, 
using  a  formalism  developed  by  Fano,’  as  an  inter, 
ference  between  the  scattering  by  a  discrete,  one- 
phonon  state  and  a  continuum  of  states.  Such  line 
shapes  have  been  reported  previously  in  BaTi03 
by  Rousseau  and  Porto"  and  in  heavily  doped  sili¬ 
con  by  Cerdeira  et  al.'°  The  application  of  Fano’s 
formalism  to  Raman  scattering  has  been  treated 
in  detail  by  Scott.11  If  the  Hamiltonian  H  couples 
a  discrete  state  <f  with  a  continuum  ipE  according 
to 

<#*l*k>**r*.  (1) 

(fit.  \n\$g)  =  E6(E.E'), 

then  the  observed  scattering  cross  section  a  is 
related  to  the  cross  section  aE  which  would  be  pro¬ 
duced  by  the  continuum  in  the  absence  of  any  cou¬ 
pling  by 

a=|(g+t)z/(l  +  eii)|ot.  (2) 

Here  €  =  (£-£„-  5£)/r  is  a  reduced  energy  vari¬ 
able  involving  the  half-width  r  =  ir  |  |3  and  energy 

shift  5£  =  ir'lP  J  tlE'  T^(E  -  E')  (P  indicates  “prin¬ 
cipal  part  of")  of  the  “resonant”  state  resulting 
from  the  discrete -continuum  interaction.  The 
parameter  <7  =  (jrFB)'1(4>|a',Ji)/(i/<E|au|i>  depends 
on  the  matrix  elements  of  the  Raman  tensor  a(/ 
for  transitions  from  an  initial  state  |t)  to  the  con¬ 
tinuum  il>B  and  to  the  state  4>  =  <p  +  pf  tiE'  VFi (iP/ 
(£-£'),  which  is  the  discrete  state  <p  modified  by 
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no.  6.  Comparison  of  tho  experimental  X(ZZ)X 
spectrum  (solid  curve)  with  the  prediction  of  Kq.  (2) 
(dots). 


an  admixture  of  continuum  states.  If  q  and  T  are 
independent  of  E  in  a  sufficiently  large  interval 
around  £„,  the  spectrum  can  be  fitted  to  Eq.  (2) 
with  q  and  T  treated  as  (energy-independent)  fit¬ 
ting  parameters.  (Note  that  F  and  q  cannot  be 
completely  independent  of  energy  unless  6E  =  0.) 
Figure  6  compares  the  experimental  X(ZZ)X 
spectrum  near  430  cm'1  with  the  prediction  of  Eq. 
(2)  to  which  a  slowly  varying  background  has  been 
added.  A  good  fit  is  obtained  for  q  =  0.75,  T  =  3.6 
cm'1,  i  6E  =  432.7  cm'1.  The  value  of  E»  ob¬ 
tained  from  the  A’(  KF)A’  spectrum  |Fig.  3(b)  |  is 
430  cm*1,  yielding  a  shift  OK -2.1  cm*1  due  to  the 
coupling.  The  exact  nature  of  this  coupling  re¬ 
mains  to  be  determined.  If,  as  we  have  proposed, 
the  continuum  is  produced  by  disorder- induced 
first-order  scattering,  it  may  be  coupled  to  zone- 
center  modes  via  anharmonic  terms  in  //  as  dis¬ 


cussed  by  Scott"  or  via  additional  terms  which 
depend  directly  on  the  disorder. 

CONCLUSION 

The  experiments  described  here  demonstrate 
that  a  major  part  of  the  dynamical  behavior  of 
orthorhombic  KNIO,  is  one  dimensional  in  char¬ 
acter  and  consistent  with  the  formation  of  linear 
chains  of  correlated  displacements  of  niobium 
ions  as  proposed  in  connection  with  diffuse  x-ray 
scattering  results1  and  Raman  scattering  studies 
of  the  cubic  and  tetragonal  phases/’  This  further 
evidence  of  disorder  in  the  orthorhombic  phase, 
together  with  the  observation  that  the  lowest 
f?2(  TO)  mode  does  not  display  any  significant 
softening  with  decreasing  temperature,  strongly 
suggests  that  the  phase  transitions  of  KNbCL,  (in¬ 
cluding  the  ferroelectric -paraelectric  one)  have 
some  order-disorder  character.  Such  a  sugges¬ 
tion  is  not  inconsistent  with  the  fact  that  the  transi 
tions  are  thermodynamically  of  first  order,  since 
recent  studies1-'  of  systems  of  dipoles  interacting 
via  a  generalized  molecular  field  have  shown  that 
first-order  transitions  are  possible  in  such  a 
model  and  that  an  accurate  quantitative  descrip¬ 
tion  of  the  cubic-tetragonal  transition  in  BaTiO, 
can  be  obtained. 
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APPENDIX  B 


Abstract  Submitted 
for  the  Washington  Meeting  of  the 
American  Physical  Society 
March  1978 


Physical  Review 
Analytic  Subject  Index 
Member  77 . SO 


Bulletin  Subject  Heading  • 

In  which  Paper  should  be  placed 
Ferroelectrlcs 


Temperature  Dependence  of  the  Devonshire  Coefficients 
of  Ferroelectrlcs."  Y.H.  TSVO,  M.I.  BELL,  and  P. X.  RACCAHt, 
Yeshlva  P.  -  The  nonlinear  dielectric  constants  of  BdTlO- 
and  TGS  at  temperatures  above  T  have  been  amaaured  by 
means  of  a  large  dc  bias  field  and  a  small  audio  frequency 
signal.  The  coefficient  of  the  term  In  Devorshlre's  ex¬ 
press  Ion  for  the  free  energy  which  la  proportional  to  the 
fourth  power  of  the  polarization  Is  found  to  have  a  strong, 
linear  temperature  dependence  In  both  cases.  The  coeffici¬ 
ent  of  the  sixth-power  term  Is  also  a  linear  function  of 
temperature  In  TGS.  These  results  are  in  qualitative  agree¬ 
ment  with  a  recently  proposed  generalized  molecular  field 
theory  of  ferroeleetrlclty,  although  certain  anomalies  which 
appear  In  the  data  for  TGS  will  require  further  study. 

*  Supported  by  ARPA  contract  DAAB-07-76-C-1342. 
t  Present  address:  U.  of  III.,  Chicago  Circle. 
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APPEIOIX  C 


* 

THEORY  OF  PYROELECTRIC  MATERIALS 
by 

M.  I.  Bell 


Abe  tract 


Elementary  thermodynamics  and  simple  microscopic  models  are  used 

to  examine  the  properties  of  pyroelectric  materials  which  determine 

their  performance  as  infrared  detectors  and  vidicon  targets.  It  is 

shown  that  for  small-area  detectors  the  signal-to-nolse  ratio  is 

/ 

subject  to  a  theoretical  limit  dictated  entirely  by  thermodynamic 
considerations.  In  a  large  class  of  simple  microscopic  models  for 
intrinsic  (proper)  fer roelectrics ,  the  ratio  of  the  pyroelectric 
coefficient  and  dielectric  constant  (which  is  a  figure  of  merit  for 
vidicon  target  performance)  is  found  to  exhibit  an  interesting 
universal  behavior.  This  behavior  can  be  exploited  to  yield  estimates 
of  the  figure  of  merit  which  suggest  that  improvements  of  more  than  a 
factor  of  three  are  most  unlikely  using  intrinsic  ferroelectrics .  An 
examination  of  extrinsic  (improper)  ferroelectrics  using  phenomenologi¬ 
cal  models  for  the  free  energy  shows  that  the  vidicon  figure  of  merit 
is  not  subject  to  the  limitations  found  for  intrinsic  ferroelectrics, 
so  that  these  materials  offer  the  possibility  of  substantially  improved 
vidicon  performance. 


*  Presented  at  the  July,  1976  meeting  of  the  Defense  Advanced  Research 
Projects  Agency  Materials  Research  Council. 
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